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摘 要 : 为 明确 利用 无 人 机 影像 监测 玉米 群体 株 高 的 精度 及 其 影响 因素 ， 本 研究 基于 无 人 机 搭载 光学 成 像 
设备 构建 大 田 玉米 群体 数字 高 程 模型 (Digital Elevation Model, DEM), ， 研 究 不 同 生育 时 期 下 玉米 群体 株 高 监 
测 的 精度 差异 。 和 针对 3 个 玉米 品种 、8 个 播 期 处 理 构建 的 株 高 差异 化 群体 ， 用 多 旋翼 无 人 机 搭载 高 清 RGB 相 
机 和 多 光谱 成 像 设 备 ， 采 集 试验 区 高 清 RGB 和 多 光谱 影像 ， 获 取 玉 米 群 体 数字 高 程 信息 DEM 和 各 处 理 区 植 
株 高 度 ， 分 析 不 同 品 种 和 播 期 处 理 下 基于 无 人 机 和 人 工 测量 株 高 间 的 相关 关系 。 试 验 结果 表明 ， 高 清 RGB 
相机 和 多 光谱 成 像 设 备 获 取 的 DEM 均 能 反映 玉米 群体 的 高 度 差 异 。 高 清 RGB 相机 的 株 高 监测 精度 优 于 多 光 
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谱 成 像 设 备 ， 但 株 高 监测 精度 不 足 ， 难 以 反映 玉米 群体 的 较 小 株 高 差异 。 不 同 生 育 阶段 对 玉米 株 高 监测 精度 
具有 较 大 影响 ， 生 育 前 期 冠 层 尚未 全 部 覆盖 地 表 或 生育 后 期 植株 衰老 叶片 枯黄 下 垂 时 ， 受 裸露 地 表 影 响 ， 群 
体 株 高 被 严重 低估 。 本 研究 分 析 了 影响 无 人 机 搭载 成 像 设 备 监测 玉米 株 高 精度 的 因素 ， 可 为 该 方法 应 用 于 大 


田 生产 提供 借鉴 意义 。 
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1 引言 


株 高 是 表征 作物 长 势 的 重要 指标 o, EEK 
状况 监测 、 水 肥 调 控 决 策 以 及 预测 倒伏 发 生 等 方 
面具 有 广泛 的 应 用 价值 ”' 。 传 统 株 高 测定 方法 
是 以 塔 尺 等 长 度 测量 工具 对 田间 或 取样 的 整 株 样 
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本 进行 人 工 测定 一 ， 不 仅 操作 效率 低 ， 还 会 因 操 
作 规 范 不 统一 而 产生 主观 误差 。 

为 提高 株 高 测定 的 效率 和 精度 ， 许 多 研究 尝 
试 利用 机 融 视 觉 、 激 光 测 距 (Laser Triangula- 
tion) "等 方法 构建 控制 环境 下 植株 个 体高 度 
的 精准 测定 方法 ， 提 高 了 株 高 测定 的 客观 性 和 准 
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确 性 。 但 这 些 方法 需要 构建 标准 化 的 检测 环境 或 
无 法 在 田间 开放 式 环 境 下 应 用 ， 实 用 性 和 测定 效 
率 仍 存在 不 足 US, 

近年 来 ， 无 人 机 搭载 数字 成 像 设备 ” 、 
高 /多 光谱 成 像 设 备 ""、LiDAR |"! 等 激光 设备 
构成 了 灵活 、 高 效 的 农田 监测 系统 ， 已 发 展 成 为 
作物 群体 生长 及 营养 监测 的 重要 手段 。 在 植株 高 
度 测定 方面 ，LiDAR 设 备 获取 的 点 云 密度 高 、 精 
度 好 ， 具 有 明显 的 优越 性 "2 。 但 LiDAR 设 备 单 
价 高 、 自 重大 、 对 搭载 的 无 人 机 要 求 高 、 数 据 处 
理 难度 大 等 问题 制约 了 其 从 科研 测定 手段 向 农业 
生产 应 用 的 拓展 。 而 高 清 RGB 相机 或 光谱 成 像 
仪 等 成 像 设备 与 无 人 机 结合 构建 的 低空 监测 平 
台 ， 可 利用 运动 重建 (Structure from Motion, 
SfM) 算法 获取 被 摄 地 物 的 稀 玻 点 云 ， 继 而 重 构 
三 维 模 型 "”， 最 终 获 取 被 摄 地 物 的 数字 表面 模 
型 (Digital Surface Model, DSM), ， 并 应 用 于 多 
种 作物 株 高 测定 。 为 消除 地 形 因 素 对 株 高 测定 结 
果 的 影响 ， 需 利用 播种 前 无 作物 的 土壤 DSM 与 
株 高 测定 时 期 的 作物 群体 DSM 进行 差 值 计 
算 2 ”。 此 外 ， 成 像 设备 的 影像 畸变 以 及 SfM 算 

法 影响 作物 群体 三 维 重 构 和 DSM 的 解 算 精 度 , 
利用 地 面 控制 点 及 其 高 精度 地 理 坐 标 可 对 影像 进 
行 几何 校正 ， 从 而 提高 DSM 的 精度 TU. 

作为 一 种 稀 植 作物 ， 玉 米 植株 生长 发 育 和 成 
熟 衰老 的 形态 结构 改变 会 引起 其 群体 及 冠 层 表面 
结构 的 显著 变化 ， 对 无 人 机 成 像 设备 重 构 作 物 群 
体 三 维 结构 以 及 解 算 DSM 造 成 困难 。 本 研究 通 
过 无 人 机 搭载 RGB 和 光谱 成 像 设 备 监测 玉米 群 
体 株 高 ， 分 析 不 同 玉米 品种 、 不 同 生育 阶段 以 及 
成 像 设备 分 辩 率 等 因素 对 株 高 监测 精度 的 影响 ， 
评价 在 大 田 玉米 生产 中 利用 无 人 机 影像 进行 株 高 
差异 性 监测 的 可 行 性 ， 为 研究 构建 适 于 大 田 作物 
实际 应 用 的 快速 、 无 损 且 灵活 性 高 、 低 成 本 的 株 
高 监测 方法 提供 参考 借鉴 。 


2 材料 与 方法 


2.1 试验 设计 

试验 于 2018 年 在 中 国 农业 科学 院 河南 新 乡 
综合 试验 站 (35° 07'52"N, 113?4557"E) 内 进 
行 。 试 验 地 块 平整 ， 海 拔高 度 在 61.84~61.87 m 
之 间 (图 1，2019 年 5 月 11 日 获取 的 试验 区 裸 土 
DEM). MÆ 139 (FK139) 、 京 农 科 728 
(JNK728) 和 郑 单 958 (ZD958) 三 种 株 高 不 同 
玉米 品种 ， 设 置 2018 年 内 4 月 20 日 (B1)、4 月 
30H (B2),5H10H (B3), 5H24H (B4), 
6H3H (B5), 6H13H (B6, 6H23H (B7) 
和 7 月 3 日 (B8) 8 个 播 期 ， 共 计 24 个 处 理 ， 大 
区 种 植 ， 每 个 小 区 面积 162 m (22.5 m x 7.2 m)， 
行距 0.6 m， 种 植 密度 75,000 株 /hm*。 田 间 管 理 
同 大 田 。 品 种 及 播 期 处 理 布置 见 图 2。 


| | 5i 
| m 66.26 /m 
- 61.29 /m 
图 1 2019 年 5 月 11 日 获取 的 试验 地 块 (35°%07'52"N， 
113°45'37"E) 数字 高 程 模型 DEM 


0 10 20 40m 
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Fig. 1 DEM of the experimental field (35?07'52"N,, 
113?45'37"E) obtained on May 11, 2019 


2.2 基于 无 人 机 图 像 的 群体 株 高 获取 


2.2.1 无 人 机 图 像 获取 
在 晴朗 无 云 、 风 力 小 于 2 级 的 天 气 条 件 下 ， 
于 2018 年 7 月 25 日 和 8 月 27 日 上 午 10:00 一 13:00， 
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图 2 3 个 玉米 品种 及 8 个 播 期 处 理 的 田间 布置 

Fig. 2 Field layout of3 maize varieties and 8 Sowing date 
treatment 

利用 DJI M600 Pro 六 旋 辟 无 人 机 搭载 高 清 数码 相 
机 (SONY o7I， 有 效 像素 2430 万 ， 最 高 分 辩 率 
为 6000 px x 4000 px, SONY Sonnar T* FE 35mm 
£/2.8 ZA ERER) 和 多 光谱 相机 (Micasense 
Rededge MX， 最 高 分 辩 率 为 1280 px x 960 px), 
飞行 高 度 设 定 为 70m， 航 向 重 琶 率 85%， 航 间 重 
羞 率 80%， 获 取 24 个 小 区 的 高 清 数码 RGB 影像 
和 多 光谱 影像 。 其 中 ,7 月 25 日 各 处 理 处 于 拔节 
期 和 灌浆 期 之 间 ， 此 时 各 播 期 处 理 群 体 株 高 差异 
最 大 ; 8 月 27 日 所 有 处 理 均 已 吐 丝 ， 各 播 期 处 理 
群体 株 高 差异 最 小 。 
2.2.2 ”无 人 机 图 像 处 理 

用 Agisoft PhotoScan Professional 软件 对 无 人 
机 图 像 进行 拼接 处 理 ， 生 成 试验 田 的 RGB 和 多 
光谱 数字 高 程 模型 (Digital Elevation Model, 
DEM). RGB 图 像 依据 DJIM600Pro 无 人 机 飞行 
控制 系统 记录 的 拍摄 时 定位 定 姿 系统 (Position 
and Orientation System, POS) 数据 进行 影像 拼 
接 ; 多 光谱 数据 则 依据 Rededge-M 设 备 配 置 的 全 


球 定位 系统 (Global Positioning System, GPS) 
导航 数据 进行 影像 拼接 。 

影像 拼接 流程 为 : 首先 根据 空间 姿态 数据 和 
相 邻 影像 的 特征 像素 对 图 像 像 素 进行 配对 ， 并 进 
行 姿 态 优化 ; 然后 逐步 生成 地 面 影像 的 稀 疏 点 云 
和 具有 精确 空间 信息 属性 的 密集 点 云 ; 再 构建 监 
测 区 域 表面 的 几何 结构 和 DSM， 解 算 为 不 同 处 
理 的 玉米 群体 表面 海拔 高 度 信息 。 各 处 理 区 的 植 
株 高 度 以 该 小 区 DSM 提取 高 程 点 的 集合 ， 去 掉 
15% 的 高 值 与 15% 的 低 值 数据 点 ， 以 剩余 70% 的 
数据 点 平均 高 程 值 作为 该 小 区 的 植株 高 度 提取 数 
据 。 利 用 ArcMap 软件 处 理 DSM 数据 并 做 图 。 


2.3 各 小 区 群体 株 高 人 工 测 量 方法 


在 不 同 处 理 的 试验 区 内 随机 选取 20 株 玉米 
作为 定 株 测定 样本 ， 均 匀 分 布 于 处 理 小 区 内 。 使 
用 伸缩 尺 分 别 在 获取 无 人 机 图 像 后 的 当天 测量 株 
高 。 未 到 吐 丝 期 的 播 期 处 理 ， 测 量 植株 自然 状态 
下 最 高 点 到 地 表 的 距离 为 株 高 ; 已 到 吐 丝 期 的 播 
期 处 理 , 测 量 植株 顶端 到 地 表 的 距离 为 株 高 。 以 
20 株 样本 的 平均 值 近似 作为 该 处 理 群 体 的 株 高 。 
各 处 理 在 测量 株 高 时 所 处 的 生育 时 期 如 表 1 
所 示 。 

7 月 25 日 ，B1~B5 处 理 下 ， 各 品种 均 达 到 
吐 丝 期 ， 其 中 Bl 处 理 早熟 品种 FK139 吐 丝 后 41 


表 1 获取 株 高 数据 时 各 播 期 所 处 生育 时 期 
Table 1 The growth period of each sowing date when the 


plant height was obtained 


7 月 25 日 8 月 27 日 
FK139 JNK728 ZD958 FK139 JNK728 ZD958 


播 期 


Bl AS41d AS36d AS34d AS74d AS69d AS67d 
B2 AS3Id AS29d AS26d AS64d AS62d ASS9d 
B3 AS21d AS18d AS16d AS54d ASSId AS49d 
B4 AS10d AS7d AS5d AS43d AS40d AS38d 
B5 AS5d ASId Silking AS38d AS34d AS33d 
B6 V14 V12 V12 AS27d AS24d AS22d 
B7 V7 V8 V7 AS17d ASI2d ASIOd 
B8 V6 V6 V6 AS8d AS4d AS2d 


注 :“AS*d” 表 示 吐 丝 后 * 天 ,“V*” 表 示 * 展 叶 期 
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天 ， 部 分 植株 叶片 枯萎 下 垂 。 而 B6~B8 处 理 ， 
三 个 品种 仍 处 于 不 同 的 营养 生长 阶段 ， 株 高 差异 
较 大 。8 月 27 日 ， 全 部 播 期 处 理 下 各 品种 均 达 吐 
丝 期 ， 其 中 早熟 FK139 在 B1~B4 人 处 理 下 多 数 植 
株 叶 片 枯萎 下 垂 。 持 绿 性 较 好 的 品种 ZD958 在 
B1 处 理 下 植株 叶片 枯萎 下 垂 。 


2.4 分 析 与 统计 方法 


单 次 影像 生成 的 DEM 数 据 是 玉米 冠 层 顶部 
的 海拔 信息 ， 包 括 植株 高 度 和 地 面 海拔 高 度 两 部 
分 。 由 于 试验 地 块 平坦 无 明显 起 伏 ， 各 处 理 地 面 
海 拨 高 度 基 本 一 致 ， 因 此 处 理 间 冠 顶 海拔 高 度 差 
即 为 株 高 差 。 但 无 人 机 航 迹 中 记录 的 海拔 信息 为 
气压 计 高 度 ， 造 成 各 次 采集 时 海拔 高 度 不 一 致 。 
为 使 无 人 机 测量 的 冠 层 高 度 与 人 工 测 量 株 高 度 具 
有 可 比 性 ， 将 基于 无 人 机 的 冠 层 顶 部 高 度 与 人 工 
测量 的 群体 株 高 数据 均 按 公式 (1) 进行 归 一 化 
处 理 : 

H - minH 


Hn = i (1) 
maxH - minH 


其 中 ，Hn 为 归 一 化 后 的 基于 图 像 的 冠 层 高 
度 或 人 工 测 量 群体 株 高 ， 为 无 量 纲 参数 ; HON 
进行 归 一 化 处 理 的 数据 m; mind A ll E AY dc 
小 株 高 ，m; max 万 为 测定 的 最 大 株 高 ，m。 无 人 
机 监测 冠 顶 相对 高 度 记 为 Hn， 人 工 测量 群体 相 
对 高 度 记 为 Hn,。 归 一 化 处 理 后 ， 各 处 理 的 冠 层 
株 高 均 为 相对 高 度 ， 数 据 区 间 范 围 为 [0，1]。 
H SPSS 20.0 和 Microsoft Excel 2010 软件 进行 数 
据 ANOVA 方 差分 析 和 作 图 。 


3 结果 与 分 析 


3.1 不 同 处 理 下 群体 株 高 差异 及 图 像 表现 


在 不 同 玉米 品种 和 播 期 处 理 ， 各 小 区 人 工 实 
测 平均 株 高 结果 如 表 2 所 示 。 各 处 理 群体 株 高 整 


表现 出 Bl 和 B2 的 群体 株 高 明显 低 于 B3~B5; 
而 B6~B8 人 处 理 分 别处 于 V6~V14 的 营养 生长 阶 
段 ， 其 株 高 低 于 Bl1~B5 处 理 ; B6~B8 处 理 间 
群体 株 高 差异 大 ， 且 品种 间 的 株 高 差异 小 于 播 
期 间 的 差异 。 而 在 8 月 27 日 ， 各 处 理 均 达到 吐 
丝 期 ( 表 1)， 也 能 观测 到 不 同 处 理 间 冠 层 株 高 
存在 差异 。 

表 2 7 月 25 日 和 8 月 27 日 两 次 监测 时 人 工 测定 的 各 小 区 

平均 株 高 


Table 2 Average plant height of two treatments measured 


artificially on 25 July and 27 August 
单位 : m 
7 月 25 日 8 月 27 日 
FK139 JNK728 ZD958 FK139 JNK728 ZD958 


播 期 


Bl 2.3b 2.5b PAK 2.3d 2.5d 2.1d 
B2 2.2c 2.6a 2.26 2.le 2.5d 2.2cd 
B3 2.4a 2.6a 2.5b 2.4bc 2.6cd 2.5b 
B4 2.5a 2.7a 2.4b 2.5ab 2.7bc 2.5b 
B5 2.5a 2.8a 2.6a 2.5a 2.9a 2.6a 
B6 1.9d 2.0c 1.7d 2.4bc 2.8ab 2.5b 
B7 1.2e 1.0d 0.9e 2.4c 2.6cd 2.2c 
B8 0.3f 0.3e 0.3f PAS 22e 1.9e 
注 : 同一 列 内 不 同 字母 代表 0.05 的 差异 显著 水 平 
将 原始 DSM 图 像 经 过 小 区 手 膜 后 得 到 各 处 
理 小 区 的 DSM KR (图 3)， 获 得 的 两 期 RGB 
Ef (图 3 (a) 和 图 3 (c) ) 和 多 光谱 图 像 
(图 3 (b) 和 图 3 (d) ) 中 各 处 理 群 体高 度 的 变 
化 均 与 人 工 测量 的 各 处 理 株 高 变化 ( 表 2) 具有 
相同 趋势 ， 表 明基 于 无 人 机 图 像 方法 可 反映 群体 


高 度 差异 变化 。 
3.2 无 人 机 影像 监测 群体 高 度 与 人 工 实 测 
株 高 间 的 关系 
将 RGB 图 像 和 多 光谱 图 像 中 各 人 处理 平均 高 
度 与 实测 株 高 进行 归 一 化 ， 对 各 处 理 基 于 无 人 机 


体 呈 现 出 从 北向 南 先 增高 后 降低 的 变化 趋势 。 同 
一 品种 不 同 播 期 处 理 间 因 所 处 生长 期 不 同 ， 其 群 
体 株 高 有 显著 差异 。 在 7 月 25 日 ， 三 个 玉米 品种 
B1~B5 已 抽雄 吐 丝 ， 株 高 已 固定 ， 且 各 品种 均 


图 像 的 归 一 化 平均 高 度 Haw 与 人 工 测 量 株 高 
归 一 化 Hn 进行 相关 分 析 。 结 果 如 图 4 (134 
页 ) 所 示 。 由 于 在 7 月 25 日 各 播 期 处 理 所 处 生 
育 时 期 不 同 ， 各 处 理 群体 高 度 之 间 差 异 较 大 ， 


ChinaXiv 合 作 期 刊 


133 


202302.00195v1 


chinaXiv 


mz: 
2z 


估算 植株 高 度 1m 


高 :177.8 79.1 
40m " 0 10 20 40m H 
{f& : 171.5 753 


(a)7 月 25 日 RGB 设备 获取 的 试验 区 DSM 
N 


(b)7 月 25 日 多 光谱 设备 获取 的 试验 区 DSM 
N 


EE. en 
高 :98.0 高 :78.4 
i. 392.5 i. SERI 
(c)8 A 27 H RGB 设备 获取 的 试验 区 DSM (d)8 月 27 日 多 光谱 设备 获取 的 试验 区 DSM 
图 3 两 次 测定 时 间 RGB 和 多 光谱 设备 获取 的 试验 区 DSM 


Fig. 3 The DSMs of experimental field constructed using RGB and multispectral equipment 

during two monitoring processes 
在 图 像 中 这 种 差异 表现 明显 ，RGB 图 像 平均 高 
度 归 一 值 Hw 同人 工 实测 株 高 归 一 值 Hn 在 7 
H25H (图 4 (a) ) 极 显 著 相 关 (r= 0.839) ; 
在 8 月 27 日 ， 各 播 期 处 理 均 达 到 其 株 高 的 最 大 
值 处 理 间 群 体 株 高 差 变 小 ， 基 于 RGB 影像 的 冠 


层 高 度 Any yy SA TWE AY BRIS Hue IR] TH 
系 仍 达 显著 水 平 (图 4 (c) ), 但 相关 系数 较 
低 ， 仅 为 0.461。 基 于 多 光谱 影像 的 冠 层 高 度 与 
人 工 测 量 株 高 间 也 存在 类 似 关 系 ， 其 中 ，7 月 25 
日 两 者 达到 极 显 著 相 关 (图 4 (b) )， 但 相关 系 
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数 (x =0.587) 低 于 RGB 图像， 说 明 除 处 理 差 异 
|. fü BE RITE s] Pop DU ets HE; 而 8 月 
27 日 ， 各 处 理 株 高 达到 最 大 值 ， 处 理 间 的 株 高 差 
异 较 小 ， 同 时 又 受到 传感器 精度 的 影响 ， 基 于 多 
光谱 图 像 的 高 度 监测 结果 表现 较 差 ， 与 人 工 测量 
高 度 间 的 相关 性 未 达 显著 水 平 (图 4 (d) )。 说 


A FK139 + JNK728 
1.0 ow 
EA 
0.8 Qe 
"T 
e oe 
06 ow Ok 
s 0. ve 
3 an) 
= Ee 
= 04 ze 


* P 
Sí t FK139 R? = 0.6699 
JNK728 R? — 0.9098 
ZD958 R? — 0.8938 


0.6 0.8 1.0 


(a)7 A 25 8 RGB 相机 监测 与 实测 株 高 相对 值 的 关系 


1.0 7 p ea 9 
ie IA 
os f -9 
oc Tet 
0.6 En 
A t ee o” 
$ xus o-9 
= A 
Zoo4€ 
P FK139 R? = 0.2287 
0.2 be JNK728 R = 0.5584 
ZD958 R? = 0.636 


0.0 0.2 0.4 0.6 0.8 1.0 


(c)8 A 27 A RGB 相机 监测 与 实测 株 高 相对 值 的 关系 


明 在 不 同 处 理 间 株 高 差异 较 大 的 情况 下 ， 利 用 高 
清 数码 相机 (RGB) 和 多 光谱 相机 可 实现 群体 高 
EWE, HRGB 图 像 测 量 精度 优 于 多 光谱 相机 。 
但 当 处 理 间 高 度 差异 变 小 时 ， 两 种 相机 测量 群体 
高 度 的 误差 将 增 大 ， 特 别 是 多 光谱 相机 更 为 


O ZD958 
1.0 FK139 R = 0.2358 a 
JNK728 R? = 0.5993 9 A 


ZD958 R? = 0.5435 


(b)7 月 25 日 多 光谱 相机 监测 与 实测 株 高 相对 值 的 关系 
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FK139 R? = 0.0936 
JNK728 R? — 0.0544 
ZD958 R? = 0.1723 
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(d)8 月 27 日 多 光谱 相机 监测 与 实测 株 高 相对 值 的 关系 


图 4 无 人 机 监测 植株 高 度 与 人 工 实 测 株 高 的 关系 


Fig. 4 Correlation between height of normalized image and height of artificially measured plant 


影响 图 像 监测 高 度 的 因素 分 析 


在 品种 和 播 期 处 理 下 ， 各 区 株 高 具有 显著 差 
5 ( 表 2)。 对 两 次 测量 的 各 处 理 基 于 图 像 的 冠 层 
高 度 归 一 化 值 届 rusw 和 各 处 理 人 工 测量 株 高 归 一 
化 值 Hn 的 相关 分 析 表 明 ( 表 3)，RGB 图 像 监 
测 结果 与 实际 株 高 的 相关 程度 均 高 于 多 光谱 图 
像 ， 且 两 个 测量 时 期 间 以 及 不 同 品种 间 表 现 出 差 
异 。 在 7 月 25 日 测量 的 RGB 图 像 中 ，JNK728 和 
ZD958 品种 的 Hnw 与 Hn, 间 均 为 极 显著 相关 ， 


相关 系数 7 均 达 到 0.9 以 上 ，FK139 显 著 相 关 , r 


达到 0.818。 而 多 光谱 图 像 中 群体 归 一 化 高 度 
Hnww 与 人 工 实测 株 高 Hn, lH), (X JNK728 和 
ZD958 显著 相关 (r4 3) A 0.774 和 0.737)， 
FK139 则 相关 性 不 显著， 相关 系数 7 为 0.486。 在 
8 月 27 日 的 测量 中 ，RGB RURKI Hnuw 5 Hno 
除 ZD958 达到 显著 相关 (x =0.797) 外 ， 其 余 两 
个 玉米 品种 均 不 显著 ,而 基于 多 光谱 图 像 的 
Hnww 与 Hnw 间 ， 三 个 玉米 品种 均 未 达到 显著 相 
关 。 表 明 玉 米 品种 类 型 与 其 发 育 时 期 影响 两 种 相 
机 对 群体 株 高 的 测量 ， 同 时 也 说 明 利 用 RGB 图 
像 监 测 群 体高 度 对 品种 和 发 育 时 期 的 适应 性 
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ES 
表 3 不 同 玉米 品种 归 一 化 图 像 高 度 与 人 工 实测 株 高 特征 
参数 
Table 3 Normalized image height of different cultivars and 


characteristic parameters of artificially measured plant height 


相关 系数 
测量 日 期 品种 = 
RGB 图 像 多 光谱 图 
FK139 0.818" 0.486 
JNK728 0.954" 0.774* 
7H25H 
ZD958 0.945" 0.73T 
—— 0.900" 0.672" 
FK139 0.478 0.306 
JNK728 0.549 0.272 
8H27H 
ZD958 0.797" 0.415 
一 一 0.604” 0.332 


注 : *、** 分 别 表示 在 0.05、0.01 水 平 上 差异 显著 


不 同 播 期 玉米 群体 株 高 的 差异 实质 上 是 玉米 
处 于 不 同 发 育 阶 段 或 因 生长 发 育 处 于 不 同 光 温 条 
件 造成 的 植株 高 度 差异 ， 这 种 差异 通过 图 像 方法 
可 很 好 地 表达 出 来 (图 3)。 在 7 月 25 日 获取 的 不 
同 播 期 条 件 下 RGB 图 像 结 果 与 人 工 测 量 株 高 极 
显著 相关 ， 相 关系 数 达 到 0.9。 而 在 8 月 27 日 ， 
基于 图 像 的 群体 高 度 与 人 工 测 量 株 高 间 归 一 化 相 
关 性 减弱 。 通 过 比较 两 次 测量 时 期 不 同 处 理 所 处 
的 发 育 时 期 发 现 ，7 月 25 日 测量 ， 各 处 理发 育 时 
期 处 于 V6 一 一 吐 丝 后 41 天 ， 植 株 均 未 成 熟 ， 冠 
层 叶 片 均 为 绿色 ; 8 月 27 日 测量 ， 各 处 理 处 于 吐 
丝 后 2 天 到 吐 丝 后 72 天 ， 其 中 有 8 个 处 理 处 于 吐 
丝 后 50 天 以 上 ， 即 8 个 处 理 已 成 熟 ， 而 成 熟 期 玉 
米 叶片 枯黄 较 多 ， 可 能 影响 测量 结果 ， 也 是 导致 
多 光谱 相机 测量 效果 变 差 的 主要 原因 。 说 明 测量 
时 期 对 基于 图 像 的 植株 高 度 测量 有 显著 影响 。 


4 讨论 


利用 无 人 机 搭载 光学 成 像 设 备 进行 作物 生长 
和 营养 诊断 是 当前 精准 农业 研究 的 热点 ， 植 株 高 
度 监 测 是 较 早 开展 且 技 术 较 为 成 熟 的 应 用 之 
一 上 9。SfM 是 一 种 三 维 重建 的 方法 ， 用 于 从 运 
动 (motion) 中 实现 3D 重 建 二。 利用 无 人 机 沿 


航线 运动 过 程 中 拍摄 时 间 序 列 作 物 影像 集 ， 经 过 
特征 提取 、 影 像 配 准 、 全 局 优化 以 及 数据 融合 等 
关键 步 又 ， 能 够 重 构 地 面 作物 群体 的 三 维 模型 。 
结合 无 人 机 惯性 测量 单元 (Inertial Measurement 
Unit, IMU) 和 GPS 等 机 载 系统 提供 影像 的 空间 
姿态 和 位 置信 息 ， 从 而 获得 作物 冠 层 的 
DSM! 。 相 比 其 他 株 高 监测 方法 ， 该 方法 所 
需 仪器 设备 价值 低 、 监 测 精 度 尚 可 且 无 损 高 效 ， 
有 具有 在 生产 中 大 规模 应 用 的 潜力 “1。 

本 研究 基于 当前 精准 农业 研究 中 普遍 应 用 的 
高 清 数码 相机 和 多 光谱 成 像 设 备 获取 群体 DSM 
数据 ， 比 较 不 同 品种 和 播 期 处 理 下 玉米 群体 株 高 
的 DSM 差 异 。 人 研究 结果 表明 ，DSM 信 息 反映 的 
群体 高 度 差异 与 实测 植株 高 度 差异 显著 相关 ， 由 
高 清 RGB 相机 监测 并 获得 的 群体 高 度 差异 数据 
与 实测 植株 高 度 差 的 相关 性 更 高 ， 监 测 精度 优 于 
多 光谱 成 像 设 备 。 第 一 次 (7 月 25 日) 监测 时 ， 
由 于 各 播 期 处 理 的 生育 进程 不 同 ， 群 体 的 株 高 差 
异 明 显 。 两 种 成 像 设 备 获取 的 DSM 数据 均 与 实 
测 株 高 显著 相关 。 而 在 第 二 次 (8 月 27 日 ) 监测 
时 ， 各 播 期 处 理 均 已 达 抽 雄 期 ， 植 株 高 度 不 再 变 
化 ， 此 时 群体 株 高 差异 是 由 关键 生育 阶段 气候 环 
境 因素 不 同 造成 的 ， 多 数 处 理 间 的 株 高 差异 在 
40 cm 以 内 ， 第 二 次 监测 两 种 设备 获取 的 DSM 高 
度 差 与 实测 株 高 差 之 间 相 关 性 均 不 显著 ,这 可 能 
与 监测 精度 不 足 ， 难 以 反映 较 小 株 高 差异 有 关 。 

无 人 机 监测 株 高 的 精度 受 影像 质量 、 定 位 信 
息 精 度 和 SfM 三 维 模型 重 构 算 法 等 多 方面 因素 影 
啊 。 成 像 设 备 的 图 像 分 辨 率 和 镜头 畸变 对 后 续 作 


明显 畸变 时 ， 作 物 影 像 的 细节 更 丰富 、 易 于 影像 
时 间 序 列 特征 点 匹配 。SfM 方 法 能 够 生成 更 精细 
的 点 云 信息 ， 从 而 提高 作物 三 维 模型 的 重 构 效 
果 。 在 地 面 建立 经 过 精确 测定 海拔 和 定位 坐标 的 
标志 点 体系 ， 能 够 帮助 SfM 方 法 纠正 影像 畸变 ， 
在 一 定 程度 上 消除 影像 集 POS 信息 精度 不 足 造 成 
的 DSM 高 程 数据 偏差 ， 也 是 提高 无 人 机 监测 精 
度 的 有 效 手段 “。 同 等 测定 条 件 下 ， 高 清 数码 
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相机 具备 像素 分 辩 率 更 高 的 精度 优势 ， 其 估算 的 
群体 株 高 相关 性 优 于 多 光谱 影像 分 析 估 算 的 群体 
株 高 。 基 于 DSM 提取 的 群体 高 程 包含 土壤 高 程 
和 作物 群体 株 高 两 部 分 信息 ， 在 高 程 起 伏 较 大 的 
地 块 ， 结 合 裸 土 DEM 和 群体 发 育 后 的 多 时 相 
DSM 进行 插值 分 析 是 获取 株 高 信息 的 必要 步 又 。 
但 影像 POS 信息 纪录 的 高 程 数据 是 无 人 机 获取 的 
气压 计 高 度 ， 不 同 监测 日 期 的 气压 计 高 度 差 异 较 
大 ， 如 无 地 面 标志 点 辅助 矫正 则 难以 准确 匹配 土 
IX DEM 与 作物 DSM 之 间 的 高 程 差 ， 造 成 株 高 监 
测 的 困难 室 。 本 研究 针对 平坦 地 块 尝 试 利 用 单 
次 监测 影像 提取 不 同 处 理 间 的 植株 高 度 相 对 差 
X. WX (7 月 25 日 和 8 月 27 日 ) 监测 数据 未 利 
用 地 面 标志 点 进行 几何 校正 。 图 3 结果 显示 ， 两 
种 监测 设备 在 两 次 监测 之 间 具 有 较 大 的 海拔 高 度 
差异 ， 特 别 是 RGB 高 清 相 机 的 监测 结果 ， 由 于 
直接 采用 无 人 机 POS 数据 生成 DSM 数据 ， 两 次 
监测 间 的 海拔 高 度 相 差 近 80 m。 而 多 光谱 成 像 
仪 影像 中 记录 的 是 GPS 测定 的 海拔 高 ， 稳 定性 相 
对 较 好 ， 但 依然 高 于 地 块 实际 海拔 高 度 10 m 以 
上 。 例 如 图 1 显示 经 过 RIK 配 准 的 试验 地 海拔 高 
度 为 61.8~-61.9m， 图 3 (b) 和 图 3 (d) 两 幅 显 
示 多 光谱 设备 监测 的 冠 层 DSM X 73.4~79.1 mo 

三 维 模型 的 重 构 精 度 与 获取 的 图 像 序列 间 匹 
配 产生 的 点 云 密度 密切 相关 。 而 受 影 像 分 辨 率 和 
SfM 算 法 精度 的 限制 ， 冠 层 上 部 器 官 的 细节 往往 
因 无 匹配 的 特征 点 从 而 被 忽视 ， 造 成 DSM 细节 
缺失 和 株 高 被 低估 的 现象 拉 。 本 研究 发 现 ,在 
生育 前 期 或 后 期 ， 因 个 体 较 小 或 叶片 枯萎 下 垂 造 
REAT, AH DSM 佑 算 群 体 株 高 的 精度 
显著 下 降 。 于 2018 年 5 月 11 日 监测 生成 的 图 1 反 
映 了 试验 地 块 的 平整 性 ， 此 时 已 有 两 个 播 期 处 理 
播种 并 出 昔 ， 但 因 幼 苗 矮 小 且 相 对 稀疏 ， 测 定 的 
土壤 DEM 在 播 期 1、2 与 其 他 未 播种 地 块 的 高 程 
几乎 无 差异 。 当 玉米 群体 生理 成 熟 后 ， 玉 米 叶片 
枯黄 下 垂 ， 冠 层 表 面 粗糙 度 增 加 ， 对 于 重 构 形 成 
DSM 造成 困难 ， 例 如 FK139 品 种 的 熟 期 较 早 ， 
在 第 一 次 监测 时 ， 播 期 1 处 理 的 FK139 小 区 植株 


高 度 估算 与 尚 在 V6 时 期 的 播 期 8 处 理 监测 高 程 
几乎 相同 。 因 此 ， 在 群体 发 育 较 小 尚未 形成 密集 
冠 层 的 生育 阶段 以 及 叶片 枯萎 下 垂 的 中 后 期 ， 利 
用 无 人 机 SfM 方 法 监测 株 高 时 会 存在 较 严 重 的 监 
测 精度 下 降 问题 。 


5 结论 


本 研究 分 析 了 在 大 田 生 产 中 应 用 无 人 机 搭载 
成 像 设备 监测 群体 株 高 的 精度 与 存在 的 问题 。 通 
过 航拍 影像 分 析 获 取 作 物 群 体 DSM， 不 同 品种 
和 播 期 处 理 区 平均 DSM 高 程 与 实测 株 高 具有 极 
显著 相关 性 。 结 果 表 明 ， 在 平坦 地 块 利用 单 时 相 
DSM 能 够 定性 分 析 群 体 株 高 差异 ,但 因 植 株 个 
体 较 小 或 茎 叶 枯 萎 的 稀 跑 群体 ， 株 高 监测 精度 明 
显 下 降 。 相 同 条 件 下 ， 像素 分 辨 率 更 大 的 成 像 设 
备 获取 的 图 像 估算 群体 株 高 与 实测 株 高 相关 性 
更 高 。 
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The Accuracy Differences of Using Unmanned Aerial 
Vehicle Images Monitoring Maize Plant Height at 
Different Growth Stages 


YANG Jin, MING Bo, YANG Fei, XU Honggen, LI Lulu, GAO Shang, LIU Chaowei, 
WANG Keru, LI Shaokun’ 


(Institute of Crop Sciences, Chinese Academy of Agricultural Sciences/Key Labrary of Crop Physiology and 
Ecology, Ministry of Agriculture and Rural Affairs, Beijing 100081, China) 


Abstract: The digital elevation model (DEM) of maize population in field was constructed by using optical imaging equipment 
mounted on unmanned aerial vehicle (UAV) to study the accuracy difference of maize population height monitoring at different 
growth stages. Three cultivars and eight sowing date treatments were set up to structure maize population with different plant 
heights. A multi-rotor UAV with high-definition digital camera and multispectral imaging sensor was used to take RGB images 
and multispectral images in the experiment area on July 25th and August 27th, 2018, which were the biggest and smallest differ- 
ences in plant height. The DEM data of maize population and canopy height were obtained with image pose correction, image 
mosaic, point cloud generation, and space reconstruction, et al. The canopy height and plant height were normalized, and the 
correlation between different cultivars and sowing date was analyzed based on UAV and manual plant height measurement. The 
feasibility of using DEM data of maize canopy to monitor the difference of plant height was clarified. The results showed that 
the height difference of maize population could be reflected by the digital elevation information obtained from high-definition 
RGB camera and multispectral camera. The plant height monitoring accuracy of HD RGB camera was higher than that of multi- 
spectral camera. However, the monitoring accuracy of plant height was not enough under the ready-made image equipment and 
treatment method. So, it was difficult to reflect the smaller plant height difference of maize population. Growth stage had a great 
influence on the monitoring of maize plant height. When the canopy of early growth stage has not completely covered the sur- 
face or the leaf yellow and withered in the later stage of growth. The plant height of the population affected by the exposed sur- 
face was seriously underestimated. In this study, the effects of UAV imaging equipment on monitoring maize plant height were 
analyzed. The influence factors can be used as reference for the application of the method in field production. 


Key words: accuracy difference; UAV; monitoring plant height; maize; multispectral 
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